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1.1 Introduction 

The electronic properties of the ground states of semiconductors have been studied both exper- 
imentally and theoretically for many years. Thus angle resolved photoelectron spectroscopy 
(ARPES), especially in combination with tunable synchrotron light sources, has been ex- 
tensively used to map the dispersion of bands in single crystals. However, in the past the 
experimental work has concentrated almost exclusively on the measurement of energies, and 
the theoretical valence-band structure calculations have been tested essentially only in terms 
of their predictions of eigenvalues. In contrast, relatively little attention has been given to 
the wave function of the electrons, despite the fact that wave function information provides 
a much more sensitive way of testing the theoretical model under investigation. Although 
the wave function cannot be measured directly it is closely related to the spectral momentum 
density. In the independent particle model it is simply proportional to the modulus square of 
the one-electron wave function, 

A{q,u) = ^2n jM \^^ q {q)\ 2 5{q - k - G)5{uj - e ]k ), (1.1) 

Gk 

where $(q) is the momentum-space one-electron wave function, j is the band index, k the 
crystal wave vector, and rijk and Ejk are the occupation number and energy of the correspond- 
ing one-electron state. The reciprocal lattice vector G translates the momentum q to the first 
Brillouin zone. For an interacting many-electron system the full spectral electron momentum 
density (SEMD) is given by 

A{q,u) = -G-{q^) = -- -1— (1.2) 

7T 7T [oj — h — Zj(Q, UJ)\ 
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Here G~(q,uj) is the interacting single-hole (retarded) Green's function of the many- 
electron system, £ is the self energy and h is the one-electron operator, which includes the 
kinetic energy and the Coulomb potential from the nuclei and the average of the electron 
charge cloud density (the Hartree potential). Presuming that the Green's function can be diag- 
onalized on an appropriate basis of momentum-space quasiparticle states <fij(q) (e.g. orbitals 
in atoms, Bloch waves in crystals) then for a crystal it takes the form 1 1 1 

A(q,u)= ]T \<PM)\ 2 ^(k+G)\lmG^{k,u). (1.3) 

In the absence of electron-electron interactions the non-interacting Green's function is simply 
a delta function and eq. J 1 . 3I > reduces to eq. Jl .2I >. The interacting SEMD contains much more 
information than simply the band dispersion. The main feature describes the probability of 
quasiparticle in band j having momentum k and energy oj. The center of the quasiparticle 
peak is shifted with respect to the one-electron energy Ejk and it acquires a width due to the 
finite quasiparticle lifetime. In addition electron correlation effects can give rise to significant 
satellite structures. 

The full SEMD can be measured by electron momentum spectroscopy (EMS) G]I21> in 
which the energies Eq, E% and E% and momenta ko, fci and k2 of the incident (subscript 
0) and two outgoing electrons (subscripts 1 and 2) in high-energy high-momentum-transfer 
(e, 2e) ionizing collisions are fully determined. From energy and momentum conservation 
one can determine for each (e, 2e) event the binding (or separation) energy of the ejected 
electron 



lu = E -E 1 -E 2 , (1.4) 
and the recoil momentum of the ionized specimen 

q = ki + k 2 - k . (1.5) 
The differential cross section is given by \ 2\. 

a(k ,k 1 ,k 2 ,uj) = {2ir) i k ~ 1 k 1 k 2 f ee A(q,uj). (1.6) 

Here f ee is the electron-electron scattering factor, which is constant in the non-coplanar sym- 
metric high-energy (e, 2e) kinematics used in the spectrometer at the Australian National 
University |3]|4). Thus the (e, 2e) cross section is directly proportional to the full interacting 
SEMD. Since the EMS measurements involve real momenta, the crystal momentum k not 
appearing in the expression for the cross section, EMS can measure SEMDs for amorphous 
and polycrystalline materials as well as for single crystals. 

The prototype semiconductor silicon has been used as a test-bed to investigate the influ- 
ence of electron correlations on the SEMD, A(q, u>). Many first-principles calculations have 
been carried out on bulk silicon, see e.g. refs pll6ll71l8l|9| ll0lfTT1ll2l . The majority of these 
calculations are based on the GW approximation to the interacting Green's function 1131 1141 . 
The dispersion of the bulk bands in silicon has been studied with ARPES along high symme- 
try directions (see e.g. refs 1 15l ll6llT7l V However, there has been essentially no experimental 
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data available on the shapes of the quasiparticle peaks (i.e. quasiparticle lifetimes) as a func- 
tion of momentum, and the satellite density as a function of energy and momentum. These 
properties of the SEMD arise directly from electron correlation effects and provide stringent 
tests for approximations to the many-electron problem. First-principles calculations of many 
physical quantities of interest require the interacting one-particle Green's function as input. It 
is therefore important to have reliable methods for accurately calculating and testing the real 
and imaginary parts of the self energy, and hence of A(q, u>). 

There are severe difficulties in extracting the full A(q, u) from experimental data obtained 
by other techniques. In ARPES these difficulties include knowing the specifics of the transi- 
tions involved, such as the untangling of final-state effects from the initial-state ones, and the 
strong energy and momentum dependence of the matrix elements 1 1 3 1 8 1 . ARPES is also very 
surface sensitive, which can obscure details of the bulk electronic structure (see e.g. 1171 '). In 
addition there is usually a significant background underlying the photo-electron spectrum and 
this can hide any continuous satellite contributions. Certain Compton scattering experiments 
in which the struck electron is detected in coincidence with the scattered photon, so-called 
(7, erf) experiments, can in principle map out the spectral function 1191 1201 . However, the 
energy resolution is such that it is extremely difficult to resolve even the valence contributions 
from that of the core electrons 1 2 1 1 . 

In this paper we present EMS measurements of the valence spectral momentum density 
for the prototypical and most studied semiconductor, silicon, and compare the results with 
calculations based on the independent particle approximation as well as calculations based on 
many-body approximations to the interacting one-particle Green's function. In section fl~2l we 
discuss the experimental technique. The theoretical models are outlined in section IT31 The 
results are discussed and compared with the FP-LMTO calculations and the first-principles 
many-body calculations in section fl~4l and. where appropriate, with previous ARPES data. In 
section fl~4l the role of diffraction is also discussed and in the last section a brief summary and 
conclusion is given. 

1.2 Experimental Details 

An outline of the experimental apparatus (described fully in [3 4|) and the coordinate system 
is shown in Fig. 11.11 An electron gun emits a highly collimated 25 keV electron beam, which 
enters the sample region inside a hemisphere held at +25 kV. Thus 50 keV electrons impinge 
on the target sample along the z— direction, the diameter of the beam being 0.1 mm. The 
emerging pairs of electrons with energies near 25 keV are decelerated and focussed at the 
entrance of two symmetrically mounted hemispherical electrostatic analyzers. The analyzers 
detect electrons emerging along sections of a cone defined by 8 S = 44.3°, which is chosen 
so that if all three electrons are in the same plane then there is no momentum transferred to 
the target (i.e. q — 0). In the independent particle picture this corresponds to scattering from 
a stationary electron. If the electrons are not in the same plane (i.e. fa ^ fa) then there is a 
y— component of momentum with q x = q z = 0, that is only target electrons with momentum 
q directed along the y— axis can cause a coincidence event. The electrons are detected by two- 
dimensional position-sensitive electron detectors, mounted at the exit planes of the analyzers, 
which measure simultaneously over a range of energies and y— momenta (i.e. range of angles 
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Figure 1.1: Schematics of the experimental arrangements. Incident electrons of momentum 
feo along the 2— axis eject an electron from a thin self-supporting Si crystal. The scattered and 
ejected electrons emerging along the shaded portions of the cone defined by 6 = 44.3°, are 
detected in coincidence by two energy and angle sensitive analysers. In the bottom panels the 
sample is indicated as a block with sides parallel to the (010) and (001) symmetry directions. 
Thus in (b) the spectral momentum density is measured along the (010) direction. In (c) the 
density is measured along the (110) direction, as the crystal has been rotated about the surface 
normal by 45°. In (d) the crystal has been tilted by 35.3° relative to the position in (c) so that 
the density is measured along the (111) direction. 



0). 

Two pairs of deflection plates mounted inside the high-voltage hemisphere along the sec- 
tions of the cone can be used to change the effective scattering angle by up to 1°. In this 
way |3| one can select nonzero values for the x— and/or the z— component of momentum, 
the y— component always lying in the range 0-5 a.u. (atomic units are used here, 1 a.u. = 
1.89 A -1 ). This allows one to probe the full three-dimensional momentum space. The two 
double deflectors are also used to check that the measured momenta correspond to the ex- 
pected ones p] 1221 . ensuring that there are no offsets in q x or q z as a result of any small possi- 
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ble geometrical misalignments. The orientation of the target specimen can be determined by 
observing the diffraction pattern of the transmitted electron beam on a phosphorus screen. A 
specific direction of the thin crystal target is then aligned with the y— axis of the spectrometer 
so that we measure the energy-resolved momentum density along that direction. The sample 
orientation can be changed without removing the sample from the vacuum by means of the 
manipulator mounting arrangement [3 1. 

The single crystal target has to be an extremely thin self-supporting film. The initial part 
of the target preparation followed the procedure of Utteridge et al 1231 . First a buried silicon 
oxide layer was produced by ion implantation in a crystal with (100) surface normal. A crater 
was then formed on the back of the crystal by wet chemical etching. The oxide layer serves to 
stop the etching. It is then removed by a HF dip and the sample is transferred to the vacuum. 
At this stage the thickness of the thinned part of the crystal is around 200 nm. Low energy 
(600 e V) argon sputtering is then used to further thin the sample. The thin part of the crystal is 
completely transparent. The thickness is monitored by observing the colour of the transmitted 
light from an incandescent lamp placed behind the sample. The colour changes with thickness 
due to the interference of the directly transmitted light with that reflected from the front and 
back silicon-vacuum interfaces. The thinning is stopped when the thickness reaches 20 nm 
(corresponding to a grey-greenish light). The sample is then transferred to the spectrometer 
under UHV. A thin amorphous layer could be present on the backside due to the sputtering, 
whereas the front side (facing the analyzers) is probably hydrogen terminated as a result of 
the HF dip. The base pressure in the sputtering chamber was of the order of 10~ 9 Torr and in 
the spectrometer the operating pressure was 2 x 10~ 10 Torr. 

The experimental energy and momentum resolutions were discussed in detail by Vos et 
al (TJ. The full width at half-maximum (FWHM) energy resolution is 1.0 eV, whereas the 
FWHM momentum resolutions are estimated to be (0.12, 0.10, 0.10 a.u.) for the q x , q y ,q z 
momentum components respectively. 

Even for a very thin target multiple scattering by the incoming or outgoing electrons has 
to be taken into account. These fast electrons can loose energy by inelastic collisions, such as 
plasmon excitation, or transfer momentum by elastic collisions (deflection from the nuclei). 
These scattering events move real coincident (e, 2e) events to "wrong" parts of the spectrum, 
as either the energy or momentum conservation equations (eqs. 1 1.41 and [131 are used incor- 
rectly. For polycrystalline or amorphous solids these multiple-scattering events can be mod- 
eled by Monte Carlo simulations 1 24 1 . In the case of inelastic scattering the binding energy as 
inferred from eq. I1.4l will be too high. Inelastic multiple scattering events can be deconvoluted 
from the data by measuring an energy loss spectrum for 25 keV electrons passing through the 
sample. This deconvolution is done without using any free parameters 1251 . This approach is 
used in the present case. Figure fl~2l shows the energy loss spectrum obtained with the present 
silicon crystal target. For single crystals elastic scattering from the nuclei adds up coherently 
(diffraction). The change of the incoming or outgoing momenta by diffraction changes the 
outcome of the measurement by a reciprocal lattice vector. We demonstrate later in this paper 
how different measurements can be used to disentangle the diffracted contributions from the 
primary (non-diffracted) contribution. 
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Figure 1.2: The energy loss spectrum observed in the spectrometer for the 20 nm thick Si 
crystal sample with 25 keV incident electrons. The fitted curve is used to deconvolute the data 
for inelastic scattering. 

1.3 Theory 

1.3.1 Independent Particle Approximation 

The local density approximation (LDA) of density functional theory (DFT) has long been 
established as a very useful tool for investigating ground-state properties of bulk semicon- 
ductors from first principles 1 26 27 1 . The advantage of DFT for approximate calculations 
in many-body systems is that one extracts the needed information from a one-body quantity, 
the electron density n(r). Although the one-particle eigenvalues in the theory have no formal 
justification as quasiparticle energies, in practice they turn out to be surprisingly accurate 11281 . 

We employed here the linear-muffin-tin-orbital (LMTO) method 1129 1 within the frame- 
work of DFT. The LMTO method is just one of many computational schemes derived within 
the framework of the DFT. The great practical advantage of the LMTO method is that only a 
minimal basis set of energy-independent orbitals (typically 9-16 per atom) is needed to obtain 
accurate eigenvalues (band energies). In the present study we implemented a full-potential 
version (FP-LMTO) of the method 1 30 1 . We write the one-electron wave function in a crystal 
in the tight-binding representation as the Bloch sum of the localised MT orbitals: 

= E ^ E fl A Mr R t). (1.7) 

t A 

Here k is the crystal momentum, j band index, t lattice (translation) vector and R basis vector. 
The label A defines a MT orbital centered at a given site R and it comprises the site index 
R and a set of atomic-like quantum numbers. The expansion coefficients a^ fc are found by 
solving the eigenvalue problem using the standard variational technique. 

Momentum space representation of the wave function "JJjfc is given by the Fourier trans- 
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form of the Bloch function: 



*ik(q) = J e- lq r *jk(r)dr. (1.8) 

Due to the periodic nature of the charge density the only non-zero contributions to Qjk(q) 
occur when q = k + G with G a reciprocal lattice vector. 

3>M = Y, cJ Gj(q-k-G). (1.9) 



G 



The contributions c* G fe , the Bloch wave amplitudes, are expressed through the (Fourier) 
integrals: 



c G,k 



ST 1 J e- l{k+G > r -$> ]k {r)dr (1.10) 
= n- 1 J2^A^ k+G > R f e-* {k+G > r <f> A (r)dr] 

A ^ 

Here it is assumed that the wave function fyjk is normalized in the unit cell of the volume ft. 
The limits of the three-dimensional integration indicates symbolically the whole coordinate 
space. The EMD in the occupied part of the band j is proportional to the modulus squared 
Bloch amplitudes: 

^ ' G 

where rijk is the occupation number. The EMD II. Ill is normalized to the total number of 
valence electrons per unit cell: 

2 J2 j d QPj(<l)= n e ■ (1.12) 



1.3.2 Electron Correlation Models 

The hole Green's function entering eauation ll.3l can be calculated by the many-body perturba- 
tion theory (MBPT) expansion on the Bloch wave basis II. 7> . Taking the first non-vanishing 
term in the MBPT leads to the so-called GW approximation 1311 1321 . In this acronym G 
stands for the Green's function and W denotes the screened Coulomb interaction. The GW 
approximation is known to give accurate quasiparticle energies 1 14 1. However, its description 
of satellite structures is not satisfactory. In alkali metals, for example, photoemission spectra 
show the presence of multiple plasmon satellites whereas the GW approximation yields only 
one at too large an energy. This shortcoming of the GW approximation has been resolved by 
introducing vertex corrections in the form of the cumulant expansion to the Green's function 
1331 1341 l35l . This allowed the inclusion of multiple plasmon creation. As a result the calcu- 
lated peak positions of the plasmon satellites were found to be in much better agreement with 
the experiment than those predicted by the GW scheme itself 1361 1371 l38l . 
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Formally, the cumulant expansion for the one-hole Green's function can be derived as 
follows. We choose the time representation for the Green's function, drop the band index j 
for brevity, and write it as 



G(k,t< 0) = i6(-t)e- iUkt+ch ( k > t \ 



(1.13) 



where ui% is the one-electron energy and C h (k, t) is defined to be the cumulant. Expanding 
the exponential in powers of the cumulant we get 



G(k,t) = G (k,t) 



l + C h (k,t) + ^[C h (k,t)] 2 + 



(1.14) 



where Go(k, t) = i9(—t) cxp(—ioJkt). In terms of the self-energy the Green function for 
the hole can be expanded as 

G = G + G ^Go + G ^Go^Go + ... . (1.15) 
To lowest order in screened interaction W, the cumulant is obtained by equating 

G C h = G J2Go, (1.16) 
where ^ = ^2 GW — iGqW. The first-order cumulant is therefore 

/OO pOO 
dt' J dTe luJkT ^2{k,t) . (1.17) 

This is then put back into Eq. Jl . 1 3i yielding multiple plasmon satellites. The energy-momentum 
representation of the Green's function can be restored by the time Fourier transform. 



1.4 Results and Discussions 

I. 4.1 Band Structure 

We first discuss the dispersion £jk{q), i.e. the dependence of the energy of the Bloch function 
$jfc(q) of band j on its crystal momentum k = q + G. The band with the largest binding 
energy is labeled 1, the next one 2 etc. In momentum space the Bloch function with crystal 
momentum k is non-zero only at momentum values k + G with amplitude k (see eq. 11.9k 
The band structure and momentum densities obtained by Kheifets et al. 1301 in a FP- 
LMTO calculation (based on the DFT-LDA) are shown in Fig. 11.31 the bands being labeled 
as discussed above. The bands are periodic in q space, with band 1 having a maximum in the 
binding energy at the V points. However, the only V point with significant momentum density 
in band 1 is the one corresponding to zero momentum (see lower panel in Fig. II. 3> . Thus the 
function with the lowest energy is a Bloch function with k = 0, cj j fc ~ 1, and the other 

ch. k - o. 

The cut of the first four Brillouin zones of silicon along the q z — plane is shown in Fig. 

II. 41 For a free electron solid the wave functions are plane waves and in the ground state the 
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Figure 1.3: FP-LMTO calculations 1301 of the dispersion (top panel) and momentum density 
(bottom panel) of Si for different high symmetry crystallographic directions. The total momen- 
tum density is split up into the densities of the individual bands as indicated. 



occupied states are within the Fermi sphere with radius kf i.e. \k\ < kf. The intersection of 
this sphere with the q z = plane is indicated by the dashed circle in Fig. 11.41 The lattice 
potential of silicon can be viewed as a perturbation on the free electron picture, so that the 
wave functions are Bloch functions with more than one |c^ fc | 2 > 0. The semiconductor 
silicon has 8 valence electrons per unit cell, and hence the first 4 bands are fully occupied 
with one spin up and one spin down electron per band. We will now discuss the results of the 
measurement of silicon and emphasize that the electron density for band j is at its maximum 
in Brillouin zone j. 

The sample was a thin (~ 20 nm) single silicon crystal with (001) surface normal, which 
is first aligned with the z direction (i.e. aligned with fep, see Fig. 11.11 (a)). Rotating around 
the surface normal, measurements were taken with the sample (100) direction aligned along 
the y— axis, then the (110) direction and 4 intermediate directions were aligned along the 
y— axis as shown in Fig. 11.41 In all these cases the potentials on the sets of deflector plates 
were set to ensure that the measurements passed through zero momentum (corresponding to 
r(o,o,o))- The experimentally observed density distributions, together with the results of the 
FP-LMTO calculation, are shown in Fig. 1 1.51 for the 6 directions measured through r( 0jU ,o)- 
The calculations were broadened with the experimental 1 eV energy resolution and split up 
into the 4 occupied bands. 

For the measurement of momenta directed along the (100) direction (0° in Fig. II .41 the 
theory predicts bands 1 and 2 occupied. In the first Brillouin zone band 1 is occupied, changing 
abruptly to band 2 at 0.61 a.u. (see also Fig. 11.41 and left panel of Fig. 11.31 . There is no 
band gap in the dispersion on crossing the first Brillouin zone. This is due to additional 
symmetry of the diamond lattice (see e.g. |39|). After leaving the second Brillouin zone the 
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*f = 1.22 a.u. 



Figure 1.4: The cut through the reciprocal lattice of silicon along the q z = plane with the 
first 4 Brillouin zones labeled. The Brillouin zone boundaries are labeled by the indices of 
the reciprocal lattice it bisects. The dashed circle is the Fermi sphere for a free electron solid 
with the same electron density as silicon. Different measurements through the q = point are 
indicated by the dashed lines. 

calculated density drops only gradually to zero. The measured density (left panel of Fig. II. 5i 
has the same behaviour , but also shows an additional branch at smaller binding energies, 
which merges with the main feature at 1 .2 a.u. This additional branch can also be seen in the 
calculated band structure for the case where the crystal has been rotated by 8.5°, and comes 
from band 4. From the shape of the Brillouin zone shown in Fig. 1 1.41 it is clear that the 
measurement along the (100) direction just misses Brillouin zone 4. Due to finite momentum 
resolution it is obvious that the measurement will pick up contributions from this zone, giving 
rise to the extra branch in the observed intensity. 

In the (110) symmetry direction (reached by a rotation by 45° along the surface normal) 
the first Brillouin zone crossing is two planes, (the (111) and (111) planes), making an angle 
of ±54.35° with the q z = plane. Thus the band switches from 1 to 3 at the double crossing. 
This gives rise to the classic band gap behaviour, band 1 having a minimum in binding energy 
(maximum in energy) at the Brillouin zone crossing, with its density petering out after the 
crossing. Band 3 slowly increases in intensity from zero momentum up to the first Brillouin 
zone boundary, where it has a maximum in binding energy, with increased density as one 
passes through Brillouin zone 3. The next extremum in energy, which corresponds to a mini- 
mum in binding energy, is when band 3 crosses the next set of Brillouin zone boundaries, i.e. 
on leaving Brillouin zone 3, and its intensity decreases thereafter as the momentum increases. 
The calculations and measurements are in quite good agreement with each other for these 
general features. For the intermediate angles band 1 remains dominant, band 4 is prominent, 
band 2 makes a significant contribution for directions not far from the (100) direction, and 
band 3 make small contributions close to the (110) direction. 

Also shown in Fig. ll.5l (bottom panels) is the spectral momentum density obtained along 
the (111) direction. This direction was reached by tilting the (110) aligned sample over 35.3° 
(see Fig. II . lf d~)>. Here the density is due to bands 1 and 2 with a large band gap at the 
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Figure 1.5: The LMTO calculated spectral momentum densities (total plus separate band con- 
tributions) along the directions shown in Fig. 11.41 compared with the measured results. The 
calculations have been broadened by the experimental energy resolution. The bottom panel 
shows the results along the (111) direction. 



zone crossing. Again these general features of the experiment and theory are in reasonable 
agreement. Note that the dispersion in the (111) direction could be mapped over a much 
larger momentum range than expected based on the momentum density distribution (see Fig. 
1.3l bottom panel). This is due to diffraction, as is described later. 

A more detailed comparison between the calculated and measured dispersion of the bands 
in the three high symmetry directions is shown in Fig. 11.61 Here a fitting procedure was used 
to extract the peak position of the measured density at different momenta. For comparison we 
include for the (100) and (111) directions also the ARPES data of refs. (Tfjl and El The 
agreement between the calculation and the present measurement is excellent. There is a small 
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Figure 1.6: The measured dispersion in the peak density (dots) along the (100) (left panel), 
(110) (central) and (111) (right panel) symmetry directions. The full line is the LMTO calcu- 
lation. Only those branches of the band structure are shown that are expected to have non-zero 
intensity in the EMS measurement. The open circles are ARPES data from ref. 1 16 1 and the 
diamonds, squares and crosses are ARPES data from Q7|. 



deviation near the X— point in the (110) direction, with the observed gap between bands 1 
and 3 being a little smaller than given by theory. A similar deviations for the band gaps is 
observed near the L points in the (111) direction. Generally the EMS shapes are smoother 
and better defined than the ARPES ones. 

Fig. 11.71 shows the band structures measured with the (110) direction aligned with the 
y— axis at three different values of q x , i.e. 0, 0.65 a.u., and 0.87 a.u., the latter two being along 
Brillouin zone boundaries. The offsets in the x— component of momentum are produced by 
choosing suitable settings for the deflector voltages for each set of deflectors |22|. For the 
measurement with q x — q z = 0, there is a single dispersing feature with a band gap at 
q y = 0.65 a.u. due to the crossing of the zone boundary. This momentum value corresponds 
to a K point of the band structure. Since the x— and y— directions through the V point are 
equivalent, the binding energy spectra at the two K points, (q x ,qy,qz) — (0,0.65,0) and 
(0.65, 0, 0), should be identical. This is indeed the case 1221 . 

The total band width can be obtained from the zero momentum spectrum of the measure- 
ment along the (111) direction. Here both the top and the bottom of the band are clearly 
visible. The experimental separation between the two peaks is 11.85 ± 0.2 eV. The error is 
mainly due to the uncertainty in the shape of the background under the high-binding energy 
peak. The LMTO calculation gives a value of 1 1 .93 eV. 

As an alternative way of determining the valence band width we measured the valence 
band and Si 2p core level in the same experiment. This allows us to determine precisely the 
position of the valence band bottom relative to the 2p 3 / 2 component of the core line. The fit 
result gives a value of 86.65 ±0.2 eV for their energy separation. From photoemission data 
(e.g.[40|) we know that the 2p 3 / 2 binding energy relative to the valence band maximum is 
98.74 eV. In this way we obtain a value of 12.1 ±0.2 eV for the valence band width. 
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Figure 1.7: The measured densities along the lines indicated in the Brillouin zone cut. The 
a;— axis of the spectrometer coincides with the crystal (110) direction. One line goes through 
the r point (zero momentum), one has a constant q x component of 0.65 a.u. and one has a 
constant q x component of 0.87 a.u. The latter two lines coincide with Brillouin zone edges. 



1.4.2 Diffraction Effects 

In a crystal coherent elastic scattering by the nuclear sites (dynamic diffraction) of the incident 
and/or the emitted electrons can shift the observed momentum distribution by a reciprocal 
lattice vector G I4fl 1421 l43l . In the present experiment the electron momenta are all very 
large (fco = 62.1 a.u., ki t 2 = 43.4 a.u.), so that for the dominant diffraction the vectors Gi are 
all much smaller than the respective ki. Thus the Gi must be essentially perpendicular to the 
ki vector to fulfil the diffraction condition 2fe,; • G, + G\ = 0. In the case where the surface 
normal (the (100) direction) is aligned with fco the smallest Go contributing to diffraction 
are of type (±2 ±2 0) (see ref. [41 1 for a detailed discussion). The smallest vector q that 
can be accessed is 1.22 a.u. where the density is very low (see Fig. II .31 . and it appears at 
q y = 1.22 a.u. Diffraction events of the type (0 ±4 0} produce an offset in the y— component 
of momentum by 2.44 a.u. The measured y-momentum line still passes through the origin 
so that we can still measure electrons with q = 0. Thus this diffracted beam causes a weak 
replication of the main spectrum shifted by 2.44 a.u. For the outgoing electrons the diffraction 
depends on which crystal direction is aligned with the y— axis of the spectrometer 1411 . 

In the apparatus it is possible to rotate the sample about the y— axis. This does not affect 
the direction of measurement, but it does affect the directions of the incoming and outgoing 
electrons relative to the crystal symmetry directions and hence their diffraction. In general, 
rotating away from the symmetric configuration, there will be fewer small reciprocal lattice 
vectors that can contribute to diffraction, decreasing their influence. However, such a rotation 
cannot eliminate the contributions from reciprocal lattice vectors of the type (0 ±2 0), which 
are perpendicular to the incoming and outgoing beams. 

In Fig. 11.81 we show the measured energy resolved momentum densities for two high 
symmetry directions in a grey-scale plot with the LMTO band structure in the repeated zone 
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Figure 1.8: Measured momentum densities along the (110) (left) and (100) (right) directions 
in linear grey-scale with the LMTO band structure in the repeated zone scheme superimposed. 
The top panels show the high momentum parts with an enhanced grey-scale to highlight the low 
intensity diffracted contributions. 

scheme superimposed on it. The contributions due to diffracted beams can be seen particularly 
at large momentum in the top panels with their enhanced grey-scale. At q = the two spectral 
densities should be identical, however, as can be seen from the figure, at this T point there is a 
small peak in the density close to zero binding energy in the (100) direction, which is absent 
in the (110) direction. This density is due to the (111) reciprocal lattice vector (the shortest 
reciprocal lattice vector for the diamond lattice, length 1 .06 a.u.) shifting density from rVx,i,i) 
to Tn _ ,o)- Rotation by 10° around the spectrometer y— axis removes this peak at low binding 
energies 1411 . In this way, by examining the diffracted contributions in regions where the 
spectral density should be zero, it is possible to remove the diffracted components from the 
measured spectral momentum density. This is shown explicitly in Fig. 11.91 

Diffraction effects are most obvious when the (111) direction is aligned with the spectrom- 
eter y— axis as can be seen already in the bottom panel of Fig. 11.51 The measured intensity 
distribution is given by the crosses in Fig. 11.91 which shows cuts through the data at the indi- 
cated binding energies. At the valence band maximum there are five peaks visible, at q = 0, 
±1.05 a.u., and ±2.1 a.u. The dominant peaks at 1.05 a.u. are due to momentum density at 
the top of band 2 (see Fig. 1 1 .31 . the other peaks are due to diffraction by a crystal reciprocal 
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Figure 1.9: The measured intensities (crosses) along (111) direction. The solid lines give the 
momentum densities corrected for diffraction as discussed in the text. In the right panel the 
length of the reciprocal lattice vector G = (111) is indicated as well. 

lattice vector. By assuming that the diffracted components add incoherently to the real spec- 
tral momentum density, the latter can be obtained by subtraction. The resulting momentum 
density is indicated by the solid line in Fig. 11.91 The small gap between bands 1 and 2 (see 
Fig. 11.31 shows up as a reduction in density around 8 eV. 

After subtraction of the diffracted component the peak at the top of the valence band at 
q = disappears completely, indicating that it was entirely due to diffraction with a reciprocal 
lattice vector of either G — (111) or G = ( — 1 — 1 — 1). However, in the binding energy 
range near 6 eV there are clear shoulders left after subtraction. This is to be expected as near 
the Brillouin zone boundary these plane waves are expected to contribute significantly to the 
Bloch waves. The shoulder and main line that are part of the same Bloch function are shown 
in this figure, connected by a line of the length of the reciprocal lattice vector. 

1.4.3 Many-body Effects 

Up to this point we have concentrated on those aspects of the electronic structure of silicon 
that could be understood in the independent particle picture. However, even from the semi- 
quantitative grey-scale presentations in Fig. 11.51 (and Fig. II. 8t of the measured intensities 
it is clear that they do not follow the LMTO calculations in detail. The measurements show 
maximum intensity at intermediate energies or near the top of the band, whereas the theory 
predicts maximum intensity near the bottom of the band. Also the energy-widths of the mea- 
sured density distributions at a given momentum are much larger than the theoretical ones 
even though the experimental energy resolution has been included in the calculations. This is 
due to lifetime broadening of the spectral momentum density by electron correlation effects. 
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Figure 1.10: EMS spectra along the (110) direction at q y ~ near the V point (left) and 
q y ~ 0.87 a.u. near the X point (right). An identical deconvolution procedure is followed for 
both spectra. The raw data are indicated by the dots, the data deconvoluted for inelastic multiple- 
scattering contributions by the solid line, and the squares show data that have deliberately been 
'over-deconvoluted' (see text for details). 



We will now look at some binding energy spectra at appropriate momenta in more detail in 
order to explore the role of many-body effects in the electronic structure of silicon. 

Fig. ll.lOl shows two spectra obtained along the (110) direction, one at q ~ near the T 
point and the other at q ~ 0.87 a.u. near the X point. The peak near the X point is very 
much narrower and taller than the peak near q ~ 0. The latter is also quite asymmetric. As 
well as the raw data the figure includes a curve showing the data deconvoluted, in a parameter- 
free way, for inelastic scattering using only the measured energy loss spectrum (see Fi g 1 1 . 2I > 
fli . This still leaves considerable intensity at high binding energies, extending some 20 eV 
or more above the position of the quasiparticle peak. If one tries to remove all the intensity 
at binding energies above the quasiparticle peak by further deconvolution one obtains non- 
physical negative intensities at some binding energies. 

In order to describe the data more quantitatively we performed full-scale many-body cal- 
culations using the GW and Cumulant Expansion approximations to the one-hole Green's 
function (see section [T.3.2> . The results are presented in Fig ll.l II which shows raw as well as 
the deconvoluted spectra along each of the two symmetry directions compared with the first- 
principles many-body calculations. The theories, convoluted with the experimental energy 
resolution, are normalized to the data at one point only, namely to the peak of the quasipar- 
ticle structure in the common q ~ spectra, the theories having the same total density in 
these spectra. The small low binding energy peak present in the data in the (100) direction at 
intermediate q y values and not in the calculations is due to the picking up of intensity from 
band 4 as discussed earlier. 

The GW calculation gives a peak in the satellite intensity at around 34 eV at q ~ 0. 
This is not observed in the data. This is a well-known failure of the GW approximation 
1131 . The cumulant expansion calculation gives a better fit to the data, although it too gives 
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Figure 1.11: EMS spectra, raw and deconvoluted for inelastic multiple-scattering effects, at 
selected momenta along the (100) (top) and (110) (bottom) directions. The full and dashed 
lines are the results of the cumulant expansion and GW calculations, respectively, and they are 
normalized to the quasiparticle peak at q — 0. 



smaller satellite density than the experiment. Both calculations describe the main quasiparticle 
features quite well in all the spectra. In particular they give the broadening and the asymmetric 
structure of this feature reasonably accurately at the smaller momenta. It is this broadening 
which gives rise to the reduction in the peak heights at the lower momenta. However, even 
though the calculations give significant life-time broadening, they nevertheless underestimate 
the width of the quasiparticle structure, the cumulant expansion model giving a slightly better 
fit. 

The finite quasiparticle (hole) lifetimes causes broadening of the observed features. In Fig. 
I1.12l we plot as a function of the momentum the observed width (FWHM) for spectra obtained 
along the (100) direction. Near zero momentum the binding energy is a weak function of 
momentum, and the observed width should be a good indication of the life-time broadening. 
Away from zero momentum the binding energy becomes a strong function of momentum, and 
additional broadening is observed due to the finite momentum resolution of the spectrometer 
as initial states with slightly different momenta (and hence binding energy) contribute to the 
spectra. In spite of this a decrease in width is observed with increasing momentum. The 
sharpest spectra are observed where the dispersion goes again through an extremum (now 
minimum binding energy) i.e. near the X point ((100) direction, see also fig. 11.1 Oi and the 
L point ((111) direction). Here a width of 1.35 eV is observed, a width clearly dominated by 
the energy resolution of the spectrometer. Also shown in the figure are the FWHM given by 
the many-body calculations. At zero momentum the calculated width is significantly smaller 
than the observed width (even if the calculations are broadened with the spectrometer energy 
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Figure 1.12: The FWHM of the quasiparticle peaks in the (100) direction plotted as a function 
of momentum. The measurements are indicated by the points with error bars. The thin line 
is a smooth fit to the experimentally obtained FWHM. The theoretical estimate of the peak 
width is indicated by squares (GW theory) and circles (cumulant expansion theory). The filled 
symbols refer to the calculated line width whereas the open symbols include the effect of finite 
experimental energy resolution. 

resolution). The calculated width decreases more rapidly with momentum, which is, at least 
in part, due to the fact that we did not incorporate the effects of finite energy and momentum 
resolution in the calculations. 

The density of the quasiparticle structure and the total density along the (100) direction 
are plotted in fig. 11.131 Also shown are the quasi-particle densities given by the cumulant 
expansion calculations, which somewhat underestimate the satellite density. For this figure 
the satellite density was defined as the total density above the main quasi-particle feature, 
whereas in the experiment we take the intensity extending 20 eV above the end of the quasi- 
particle peak, in this way we obtain for the experiment that the percentage of the total density 
that is contained in the satellite structure decreases as the momentum increases. For instance, 
along the (100) direction, the satellite structure has essentially disappeared by momentum 
values near 1 a.u., whereas at q = it accounts for around 40 % of the total density (Fig. 
11.13k In the theory the satellite intensity remains fairly constant (at 25 % of the total density) 
but becomes more spread out over binding energy with increasing momentum. 

1.5 Conclusions 

We have measured the spectral momentum density of the prototype semiconductor Si along 
the three high symmetry directions, and several cuts between the (100) and (110) directions, 
all going through zero momentum (the rmo.o) point). In addition several measurements were 
made along the (110) direction but with the cuts shifted from F( 0i o) point to coincide with 
nearby Brillouin zone edges. 
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Figure 1.13: The energy integrated electron momentum density along the (100) direction. The 
open squares give the total density, and the solid circles the quasi-particle density. Also included 
is the total density given by the LMTO calculation (solid line) and the corresponding quasi- 
particle intensity, as predicted by the cumulant expansion calculations (dashed line). 

The band structure, i.e. the dispersion of the position of the peak in the quasiparticle struc- 
ture, is very well described by the FP-LMTO-DFT-LDA approximation. The present EMS 
data are in better agreement with the theory than the available ARPES data. The band-width 
was determined to be 11.85 ±0.2 eV, in good agreement with ARPES data (12.5 ±0.6eV|44|) 
and the LMTO theory (1 1.93 eV). 

Diffraction of the incident and/or the emitted electron beams can shift intensity from the 
'direct' spectral density by reciprocal lattice vectors. Such effects are observed but can often 
be distinguished from the 'direct' intensity. Changing the diffraction conditions for the in- 
coming and/or outgoing electrons pinpoints which part of the measured intensity is diffraction 
related. In this way contributions due to diffraction can be removed on the assumption that 
there is no interference between the direct and diffracted electron waves. The study of any 
possible interference effects could lead to additional information on the structure of silicon 

El. 

The spectra, especially at low momenta, show very large effects due to correlations. In 
particular the quasiparticle structures are very broad due to the short hole lifetimes, and much 
of the spectral density is found in higher energy satellite structure. This satellite structure is 
quite smooth and extends to around 20 eV above the quasiparticle peak. Although the satellite 
structure dominates the density at low momentum, its density drops off as the momentum is 
increased. Similarly the lifetime broadening diminishes as the momentum increases. Compar- 
ison with first principles many-body calculations shows that the GW approximation predicts 
the main quasiparticle features reasonably well, but cannot describe the satellite structure. 
The cumulant expansion approximation also describes the quasiparticle features qualitatively 
well, but does considerably better than the GW in its description of the shape of the satellite 
density. Although it predicts the energy distribution and the momentum dependence of the 
satellite structure quite well, it nevertheless significantly underestimates the satellite intensity, 
particularly at low momenta. The fact that agreement is better but still not perfect for the more 
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detailed model is encouraging in that a quantitative comparison between the present EMS 
data and other representations of the many-body wave function should lead to new levels of 
understanding. 
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